Deriving molecular bonding from macromolecular self-assembly 
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Abstract 

Macromolecules can form regular structures on inert surfaces. We have developed a combined empirical 
and modeling approach to derive the bonding. From experimental scanning tunneling microscopy (STM) 
images of structures formed on Au(lll) by melamine, by PTCDA, and by a 2:3 mixture of the two, we 
determine the molecular bonding morphologies. Within these bonding morphologies and recognizing the 
distinction between cohesive and adhesive molecular interactions we simultaneously simulated different 
molecular structures using a lattice Monte Carlo method. Within these bonding morphologies there is a 
distinction between cohesive and adhesive molecular interactions. We have simulated different molecular 
structures using a lattice Monte Carlo method. 

PACS numbers: 68.37.Ef, 05.10.Ln 81.07.Nb 81.16.Dn 
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Self-assembly of molecules on atomically wel 
for generating two-dimensional nanostructures jl 



defined surfaces offers a bottom-up_ 
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approach 



Jj. To 



develop structures for specific purposes requires a precise knowledge of the molecular bonding, 



including molecular binding rules and the corresponding binding energies 115 



1611 . These cannot 



always be calculated ab initio because of the complex character of molecular bonding and the 



difficulty of including dispersive (van der Waals) interactions 111 711 . 



Hydrogen bond-forming molecules are particularly suitable for generating self-assembled 



structures due to the high selectivity and directionality of hydrogen bonds 119 
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2 ill , and rel- 



atively low energetics which enable equilibrium molecular configurations to be achieved at rela- 



tively low processing temperatures [|8 



molecules, with diverse binding rules, [22, 
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Using mixtures of different hydrogen bond-forming 



23J] allows the formation of a wider variety of molecular 



structures by changing molecular composition, and presents a promising approach to generating 
molecular scaffolds [24]. The relation between the properties of individual molecules and the 
characteristics of their self-assembly on surfaces remains to be understood. Our study provides a 
step towards establishing such a relation by solving the inverse problem - extracting characteristics 
of molecular interactions by analysing self-assembled structures. 

In this paper we use experimental observations of a self-assembled molecular structure to de- 
termine plausible binding rules, and then perform kinetic Monte Carlo simulations to estimate the 
binding energies. We analyse simultaneously the structural stability of the stoichiometric molec- 
ular mixture and the individual constituent molecular components. We apply this methodology to 
PTCDA and melamine molecules on Au(l 11), where the substrate has little effect on intermolec- 
ular interactions. 

Our procedure for extracting the molecular binding rules and the corresponding interaction en- 
ergies involves two major steps. In the first step we analyse experimentally observed molecular 



structures in order to extract the characteristics of the molecular arrangement for all molecular 
compositions of interest. Thereby we identify the ways molecules bind to each other, i.e. the 
molecular binding rules. In the second step we construct the model system using the binding rules 
and arranging the molecules on a grid, suitable for on-lattice kinetic Monte Carlo simulations of 
molecular structures. Thereby we extract the molecular binding energies by analysing simultane- 
ously the stability of both the single-component molecular structures and of that formed from a 
binary mixture. 

We used Au(lll) film grown on mica substrates. The samples were introduced into the ul- 
trahigh vacuum (UHV) chamber of a STM (JEOL JSTM4500S) operating at a pressure of 1(T 8 
Pa. The Au(lll) surfaces were sputtered with argon ions and annealed in UHV at temperatures 
between 600 and 800 °C typically for 30 min. PTCDA molecules were sublimated at 275 °C and 
melamine at 100 °C. Electrochemically etched tungsten tips were used to obtain constant current 
(It) images at room temperature with the bias voltage (V s ) applied to the sample. The structures of 
molecular mixture was obtained after deposition of PTCDA on Au(l 11), followed by a deposition 
of melamine with a 2:1 ratio and a post annealing at 90 °C for 10 hours. 

Molecular binding rules between identical molecules are obtained by comparing the molec- 
ular arrangements in the single-component molecular structures. The experimental images of 
the molecular structures of 3,4,9, 10-perylene-tetracarboxylic-dianhydride (PTCDA) and 1,3,5- 
triazine-2,4,6-triamine (melamine) self-assembled on Au(l 1 1) are shown in Fig. [Q FigHk shows 
the compact domains of PTCDA, which exhibit a uniform structure with a herringbone-like pat- 
tern. The 12.0 x 20.0 A 2 unit cell of the structure is rectangular and contains two molecules with 



their main axes oriented at an angle of 86° with respect to each other Il25h . Molecular bonding 



in the PTCDA structure can be characterized by two independent parameters Epp and E P p as il- 
lustrated in Fig. [Tg and lh, respectively. Melamine molecules on Au(lll), according to Fig. Id, 





FIG. 1: (a-c) PTCDA domain: a) STM-image on Au(lll) surface (14xl0nm 2 ; V s = +1.5 V, \ t = 0.4 
nA) , b) simulation result displayed on a sheared hexagonal lattice, c) molecular ordering, (d-f) Melamine 
domain: d) STM-image (14x 10nm 2 ; V s = -1.0 V, \ t = 0.5 nA), e) simulation shown on a simple hexagonal 
lattice, f) molecular ordering. The unit cell is outlined in blue. (g,h,i) Plausible molecular bonds occurring 
in structures c and f. In the molecule 3D representation, gray balls are carbon atoms, red balls are oxygen 
atoms, white balls are hydrogen atoms and blue balls are nitrogen atoms. 

form domains of chiral structure and hexagonal symmetry with the lattice parameter of 9.8 A. 



This arrangement is stabilized by a double hydrogen bond [26], as illustrated in Fig. [Qi. The 
corresponding energy E mm for the melamine-melamine bonding is found to be E mm = 0.45 eV 



021 



Details of the melamine-PTCDA molecular bonding are obtained by analysing the structure 
of the melamine-PTCDA mixture with the composition 2:1. This large scale structure is formed 
of ordered molecular stripes, as shown in Fig. |2k. Each stripe is composed of a single PTCDA 
molecular row and a double row of melamine molecules, Fig. [2Jd. The PTCDA molecular axis is 
rotated by 50 ° with respect to the stripe line thereby making this structure chiral. The unit cell 
outlined in Fig. |2b in blue has a parallelogram shape, with an angle of 85 ° with 10.0 A (the 
periodicity along the PTCDA rows) and 19.9 A parameters (PTDCA-PTCDA separation across 
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FIG. 2: Mixed PTCDA and melamine domain, (a) STM-image on Au(lll) surface (80x60nm 2 ; V s = -1.2 
V, It = 0.2 nA) , (b) Close up (14x8 nm 2 ; V s = -1.5 V, I t = 0.1 nA) , (c) simulation result shown on a sheared 
hexagonal lattice, (d) model of the molecular ordering, (e-f) plausible molecular bonds occurring in d. 

2 melamine rows). Fig. [2]3 shows the molecular arrangement as observed in Fig. |2b- Molecular 
interactions in the PTCDA-melamine structure can in general be characterized by two parameters 
Em P and Em P corresponding to physically reasonable hydrogen bonding as illustrated in Fig. [2hf. 

The molecular binding energies E PP \ E p 2 J , Em P and Em P can now be estimated by studying 
the stability of the self-assembled s ing 



e-components and the mixed molecular structure using the 
kinetic Monte Carlo methodology H27I1 with an underlying hexagonal grid defining the topology 
of molecular movements and interactions and an (N,V,T) ensemble. 

The on-lattice kinetic Monte Carlo model is an entirely topological model. It exclusively con- 
siders binding energies between molecules that are nearest neighbours in the underlying lattice, 
and does not contain information about the absolute position and orientation of molecules. There- 
fore the point of comparison between experimental and simulated structures is the topological 
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FIG. 3: Simulated PTCDA-melamine ordering depneding on PTCDA-melamine bond energies Em P = 

(2) 

E)^ v = E ad , varying from O.leV to 0.27beV. (a,b) Small bond energies between PTCDA and melamine 
with E ac i < 0.175 lead to a phase segregation. (c,d) In a narrow energy range of 0.2 < E ad < 0.225 
melamine double rows can be observed as a typical defect, (e-f) Energies E ad > 0.25 allow the formation 
of a defect-free PTCDA-melamine domain. 




FIG. 4: The double melamine row in an PTCDA-melamine domain, (a) close-up of an STM image 
(8x8nm 2 ; V s = -1.2 V, l t = 0.2 nA). (b) Scheme showing the common topology of both experiment and 
simulation. Red rectangles are PTCDA molecules and blue triangles are melamine molecules, (c) Simulated 
image using a simple hexagonal geometry. 
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correspondence (Fig. HI). The freedom to apply different geometries to a simulation while keeping 
the topology unchanged can be used for visualisation of highly ordered homogeneous molecular 
arrangements (Fig{H Fig©, where a sheared hexagonal geometry achieves a good correspon- 
dence. In other cases where the geometry is more complex due to defects or irregularities (Fig. [3j 
Fig. HJ, the simulation is presented with a simple hexagonal geometry. In each case, what matters 
is the correspondence in topology between simulation and experiment. 

In our simulations we use a binary mixture of anisotropic molecules of type 1 (trigonal vertices) 
and type 2 (linear rods) occupying one and two sites, respectively, on a 30x30 two-dimensional 
hexagonal lattice, initialized in a random configuration. Rod-like molecules exist in three different 
orientations along the symmetry axes of the hexagonal grid, while vertex-like molecules exist in 2 
different orientational configurations, due to their two-fold and three-fold molecular symmetry, re- 
spectively. Any pair of molecules that exists in one of the nearest neighbour configurations, shown 
in Fig. [[g,h,i and Fig. |2]e,f, establishes a hydrogen bond with the energies E mm , E P p 2 ' ) , Emp , 
respectively. The number of molecules (200 vertices, 100 rods) is chosen to be stoichiometric, 
complying with the ratio of the experimentally observed unit cell shown in Fig. [21 The simulation 
temperature of k B T=0.08 eV is higher than in experiment to enable short equilibration times, but 
low enough to avoid thermally generated defects. 

In our simulations we neglect the influence of the substrate on the molecular binding ener- 
gies, which is a reasonable approximation for the PTCDA and melamine molecules on Au(lll) 



[12811 . The structure of the molecular mixture can be unstable against phase separation into single- 



molecular domains. This instability is controlled by the strength of the cohesive interactions asso- 
ciated with E P p and E P p relative to the adhesive interactions associated with Emp and Em P - This 
trade-off can be simplified by assuming that E^faE^=E co h and Em P ^Emp=E adh reducing the 
number of fit parameters down to two. 



In Fig. [T]and Fig. [2] we show the predicted structures with high symmetry molecular ordering 
in order to compare with the illustrated experimental images. The pure PTCDA domain in Fig. 
Q}) and the pure melamine domain in Fig. are both the result of the binding rules and the high 
enough binding energies to ensure thermal stability of the structures. The stability of the PTCDA- 
melamine structure requires a certain range of adhesive interaction energies for given values of the 
cohesive interaction energies, as illustrated in Fig. [3] 

In Fig. [3k-f we show how the structure undergoes a change, caused by varying the PTCDA- 
melamine interaction strengths Emp=E^ p =E ad . We observe that values of E ad = 0.175 eV, E coh = 
Epp=Epp = 0.15 eV lead to the phase separation due to the weak PTCDA-melamine interactions. 
Whereas E ad = 0.2 eV creates structures with long range ordering and a sporadic but reproducible 
occurrence of double melamine row defects. Fig. |4] highlights that these double row defects in the 
PTCDA-melamine domains are observed in STM (Fig. [2^ and Fig. Hk) and are also predicted by 
our simulation (Fig. |4fc). It can be seen from Fig. HJ> that both experiment (Fig. Hk) and simulation 
(Fig. Hfc) correspond topologically. Reproducing the double-row defect in our simulations allows 
the range for the effective energy parameters to be narrowed as given in Table 1 . 

The interaction energy values we found allow all three experimentally observed structures 
shown in Fig. [Hand Fig. [2]to be simulated simultaneously, as illustrated in Fig. [Eb,e, and Fig. 
The obtained energies E ad and E coh comprise the contributions from both the hydrogen bonds and 



van der Waals interactions. The latter energies are usually in the range 0.04-0.1 eV Ill7h . which 
is around the lower limit of the variation range identified for the parameters E ad and E coh . Hence 
the major contribution to the obtained values of E ad and E coh can be associated with the hydrogen 
bond energies. 

By combining STM images of molecular ordering in a PTCDA-melamine system with on- 
lattice Monte Carlo simulations of the structural stability we have determined the molecular bind- 

8 



Type 


Bond type 


Energy 


Arrangement 


Energy range 


mm 


2xN - H---0 


E mm = 0AheV [27J 




reference value 




N — H ■ ■ ■ 


E mp = E a dh = 0.225eV 




0.2eV < E adh < 0.25eV 


m p( 2 ) 


N-H-0 


E mp = Eadh = 0.225eV 


mi 


0.2eV < Eadh < 0.25eF 


m?0) 


2 x C-H---0 


E pp = E coh = 0.15eV 




Q.leV < E^ < 0.3eV 


ppi 2 ) 


2x C-H---0 


E pp = E coh = 0.15eV 




O.leV < E coh < 0.3eV 


TABLE I: 


Summary of binding 


rules and energies in the PTCDA-melamine system. 


Abbreviation: 



mm=(melamine-melamine), pp=(PTCDA-PTCDA) and mp=(melamine-PTCDA). 

ing energies. We have decomposed the molecular interactions into adhesive and cohesive parts, 
which reduces the number of free parameters and hence the uncertainty in the energy values. The 
results provide information which can be used as a starting point for more detailed studies of 
molecular bonding using more sophisticated calculation techniques. This procedure of extract- 
ing the molecular binding rules and estimating binding energies is applicable to a wide range of 
multicomponent systems. 

Acknowledgment. The authors thank the EPSRC (EP/D048761/1 and GR/S 15808/01) for 
funding and Chris Spencer (JEOL UK) for valuable technical support. 



[1] J. V. Barth, Annu. Rev. Phys. Chem. 58, 375 (2007). 

[2] G. P. Lopinski, D. D. M. Wayner and R. A. Wolkow, Nature 406, 48 (2000). 

[3] H. L. Zhang, W. Chen, L. Chen, H. Huang, X. S. Wang, J. Yuhara and A. T. S. Wee, Small 3, 2015 
(2007). 

[4] F. Silly, A. Q. Shaw, K. Porfyrakis, G. A. D. Briggs and M. R. Castell, Appl. Phys. Lett. 91, 253109 

9 



(2007). 

[5] V. Dusastre, Nature 406, 31(2000). 

[6] P. T. Tanev and T. J. pinnavaia, Science 267, 865 (1995). 

[7] L. Chen, W. Chen, H. Huang, H. L. Zhang, J. Yuhara, and A. T. S. Wee, Adv. Mater. , in press. 
[8] J. A. Theobald, N. S. Oxtoby, M. A. Phillips, N. R. Champness and P. H. Beton, Nature 424, 1029 
(2003). 

[9] F. Silly, A. Q. Shaw, G. A. D. Briggs and M. R. Castell, Appl. Phys. Lett. 92, 023102 (2008). 
[10] D. Bonifazi, A. Kiebele, M. Stohr, F. Cheng, T. Jung, F. Diederich and H. Spillmann, Adv. Funct. 

Mater. 17, 1051 (2007). 
[11] E. Winfree, F. R. Liu, L. A. Wenzler and N. C. Seeman, Nature 394, 539 (1998). 
[12] J. V. Barth , G. Costantini and K. Kern, Nature 437, 671 (2005). 

[13] R. Otero, M. Schock, L. M. Molina, E. Lagsgaard, I. Stensgaard, B. Hammer and F. Besenbacher, 

Angew. Chem. Int. Ed. 44, 2270 (2005). 
[14] F. Silly, A. Q. Shaw, M. R. Castell and G. A. D. Briggs, ChemComm. , in press. 
[15] K. J. Franke, G. Schulze, N. Henningsen, I. Fernandez-Torrente, J. I. Pascual, S. Zarwell, K. Riick- 

Braun, M. Cobian, and N. Lorente, Phys. Rev. Lett. 100, 036807 (2008). 
[16] T. Yokoyama, S. Yokoyama, T. Kamikado, Y Okuno and S. Mashiko, Nature 413, 619 (2001). 
[17] S. Grimme, J. Comput. Chem. 25, 1463 (2004). 

[18] J. S. Swabrick, B. L. Rogers, N. R. Champness and P. H. Beton, J. Phys. Chem. B 110, 61 10 (2006). 

[19] M. M. Conn and J. Rebek, Chem. Rev. 97, 1647 (1997). 

[20] E. A. Archer, H. Gong and M. J. Krische, Tetrahedron 57, 1 139 (2001). 

[21] L. J. Prins, D. N. Reinhoudt and P. Timmerman, Angew. Chem. Int. Ed. 40, 2382 (2001). 

[22] R. E. A. Kelly and L. N. Kantorovich, Surf. Sci. 589, 139 (2005). 



10 



[23] R. E. A. Kelly, Y. J. Lee and L. N. Kantorovich, J. Phys. Chem. 109, 11933 (2005). 
[24] T. Moriuchi and T. Hirao, Chem. Soc. Rev. 33, 294 (2004). 

[25] T. Schmitz-Hubsch, T. Fritz, F. Sellam, R. Staub and K. Leo, Phys. Rev. B 55, 7972 (1997). 

[26] W. Xu, M. Dong, H. Gersen, E. Rauls, S. Vazquez-Campos, M. Crego-Calama, D. N.,Reinhoudt, I. 

Stensgaard, E. Laegsgaard, T. R. Linderoth and F. Besenbacher, Small 3, 854 (2007). 
[27] U. K. Weber, V. M. Burlakov, L. M. A. Perdigao, R. H. J. Fawcett, P. H.,Beton , N. R. Champness, J. 

H. Jefferson, G. A. D. Briggs and D. G. Pettifor, Phys. Rev. Lett., in press. 
[28] S. K. M. Henze, O. Bauer, T.-L. Lee, M. Sokolovski and F. S. Tautz, Surf. Sci. 601, 1566 (2007). 



11 



